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Abstract 

A generator with a^ order radiative correction effects, including both the initial 
state radiative corrections and the photon vacuum polarization corrections, is built 
for the full simulation of e^"e~ -^ hadrons in cross section scan experiment in the 
CN ' quarkonium energy region. The contributions from the hadron production structures 

^ • including the resonances of 1 quarkonium families and the light hadrons spectrum 

V^ , below 2 GeV as well as the QED continuum hadron spectrum are all taken into ac- 

CN I count in the event generation. It was employed successfully to determine the detection 

efficiency for selection of the events of e~^e~ -^ hadrons from the data taken in the 
^ . energy region from 3.650 GeV to 3.872 GeV covering both the ^^(3686) and V(3770) 

^^ I resonances in BES experiment. The generator reproduces the properties of hadronic 

^S^ • event production and inclusive decays well. 
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^ 1 Introduction 



'0(3770) is the charmonium resonance with the lowest mass above open charm-pair threshold. 
/\ ' It is thought to decay almost entirely to DD l\. However, it has been a puzzle for a long 
c^ ' time that some measured values of DD production cross section a^)/) at the peak of ip(3770) 
resonance fail to fit in the measured values of the cross section (7^(3770) for '?/^(3770) production 
at the peak. The discrepancy between the (7^(3770) and a^^ was historically found to be more 
than 30% |2] [SI- To understand this long-standing puzzle, it is crucial to measure the ?/'(3770) 
resonance parameters and the DD cross section at the resonance peak precisely in the same 
experiment. Experimentally, '0(3770) resonance parameters are extracted from fitting the 
observed hadronic cross sections in the energy range which covers the '0(3770) resonance in 
cross section scan experiment. However, since the ?/'(3686) resonance is closed to the -0(3770), 
and the cross section of continuum QED production exceeds the contribution of -0(3770) 
resonance the heavy overlaps between them would affect the cross section configuration in 
this energy region. In order to accurately measure the parameters of the resonances one had 
better simultaneously to deal with both the resonances and the contribution of continuum 
production based on the same data set. So, the determinations of the efficiencies for detecting 
the inclusive hadronic events point by point is necessary in the energy region covering several 
production processe. 



The determinations of efficiencies for detecting the e'^e^ — > hadrons events in such 
energy region would be more comphcate than that in the case for which only single narrow 
resonance has been concerned in the data taking, in which the initial state radiation (ISR) 
effect and the continuum contribution are negligible small, or in the case of the simple R value 
measurement in the e^e~ annihilation in the energy region far away from any resonances. 
Due to the ISR effects, the actual energy to produce final products would drop to rather wide 
region below the collision energy by photon emissions from the beam particles. The detector 
response, the composition of the different physical processes in the whole region and even the 
event selection criteria all vary with the energies and finally affect the detection efficiency at 
the collision energy points. To obtain the efficiencies at each of the different energy points 
reliably, a generator to produce the full processes of e"'"e~ — *> hadrons is developed. 

Sometimes for convenience the contribution from the vacuum polarization can be turned 
off in the event generation. The possible interferences between the continuum hadron produc- 
tion and the inclusive electro-magnetic (EM) decay mode of narrow resonances are considered 
in the generator and can be turned on/off too. 

2 Method 

2.1 Radiative correction in e+e annihilation 

Due to ISR, the actual cm. energy of the inclusive hadron system produced in the e^e~ 



annihilation is V^ = J s{l — x), where ^/s is the the collision center-of-mass energy, and 
x(l > X > 0) is a parameter relating to the total energy of the emitted photons [H]. 

Fig. n shows the differential cross section of the hadron production with respect to vs^ = 



s(l — x) for the case of that the nominal center-of-mass energy is set to 3.80 GeV. To 
a^ order radiative correction in e~^e~ annihilation, the differential cross section of hadron 
production with respect to the hard photon energy parameter x, instead of vs^ can be 
written as 

dx |i-n[s(i-x)]r ^ '^' ^^ 

where F{x, s) is a sampling function for the radiative photon energy fraction x. Neglecting 
the Z" exchange contributions in the structure function approach introduced by Kuraev and 
Fadin [5], it can be written as 

F(x, s) = f]x^-^S^+^ + (5^, (2) 

in which /3 is the electron equivalent radiator thickness, 

,3=^(i-l) 

vr 

and 

L = ln4 

ml 

where TUg is the mass of electron, a is the fine structure constant, and 

4^ TT^ 3 2' 24M 3 ' ^ ' 
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Figure 1: The differential hadron production cross section due to ISR return with respect 
to the actual cm. energy v^, da^^^{s)/d\/s' in the case that the nominal collision center- 
of-mass energy is set to y/s = 3.80 GeV. The dashed line shows the differential cross section 
for /i"^/i~ production. 
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In Eq. (^, 1/1 1 — n[s(l — x)]p is the vacuum polarization correction factor including both 
the leptonic and the hadronic terms i.e. 



n(s) = n4s) + n;(s) 



in which 

and 

where 
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and 

in which ,^ = ^^ where m^ is the mass of lepton /(/ = e, /i, r). In Eq. ((H), 

As{l - X)] = 3^ i? + E ^L,J^(1 - ^)] (6) 



is the total lowest order cross-section at y/s' = Js{l — x) in which R is the ratio of continuum 
hadron cross section to the cross section of e~^e~ — > fi~^fi~, and 



r /I M livri h^ii gg j 






is the lowest order cross section of the ith resonance, where Fg^ j, Mf , Th,i and Fj are the 
lowest order leptonic width, the lowest order mass, the hadronic width and the total width 
of the ith resonance, and 'i' denote the ith. resonance of the 1 quarkonium states and the 
states with mass below 2 GeV/c^. For wide resonances the widths Fj's and Th,i vary with 
energy depending on their decay products. 

In Fig. 121 a), the solid line shows the total lowest order cross section for inclusive hadronic 
event production in the region from 3.65 GeV to 4.0 GeV including the five processes: 
e+e" -^ Jji), e+e" -^ ^(3686), e+e" -^ V^lSTTO), e+e" -^ ?/'(4040) and e+e" ^ 7* ^ 
hadrons] while the dashed line shows the expected total cross section with the radiative 
corrections for the inclusive hadron production. These cross sections are generated with the 
generator. In Fig. |21b), the contributions of the radiatively corrected cross sections of the 
five processes are demonstrated, respectively. The energy spread for generating the events 
is set to be zero. 

The ratio of the expected observation cross section to the lowest order cross section, see 
the two lines in Fig. |21a) is the radiative correction factor for the inclusive e~^e~ -^ hadrons 
process 

'-'--^' '«' 

which is shown in Fig|21c). 

From Fig. |21 c) one can see that the radiative correction effects vary rather large around 
the resonance. Except for the ISR corrections, under the narrow resonance peak the large real 
and imaginary part of vacuum polarization correction would shift the lowest order resonance 
peak position M^ up by an amount of 

— F" -I- — f?F 

where R denotes the R value for the continuum hadron production and n'[(M°)^] denotes 
the vacuum polarization function from all the contributions of leptons and hadrons except 
for the contribution from the resonance itself at the position s = (M*^)^. And 

M = M° + AM (10) 



Sio'^r 



'10' 
10' 
10 

1 







a) 




\|/(3686) 






V(3770) 


\l/(4040) 


; Continuum QED 









lio^ 



1 



b) 



y(3686) 
1 U It Continuum QED 




J/y 




3.7 3.8 3.9 

VsGev 



10 \L 



^ 1 U ^ bare- .-.dressed Born 

a 
10 "k / 



10' 




; Continuum QED 



3.686 3.688 

VsGev 



Figure 2: The cross sections of e~^e~ anniliilation into inclusive hadron in the energy region 
from 3.65 GeV to 4.0 GeV. a) shows the total lowest order cross section (solid line) and 
the expected total cross section in the region, produced by the generator with null collision 
energy spread setting (dashed line); b) demonstrates the respective contributions of the 
expected cross sections of the five processes mentioned in the text; c) represents the radiative 
correction factor 1 + 6; and d) gives some detail information about the radiative corrections 
at the position of the sharp resonance il){3Q8Q) in the tiny energy window from 3.6845 GeV 
to 3.6895 GeV. The EM interference effect between '?/;(3686) and continuum QED is turned 
off in the event generation. 



is the observable peak position of the resonance related to the single photon exchange pro- 
cesses. One can also dress the lowest order leptonic width F^g to the physical one Fge in the 
same way, by 



^0 



:iii 



'^ |l-i?en'[(M0)2]|2' 

for calculation of da'^^^{s)/dx given in Eq. (^. The change of the total width F due to the 
vacuum polarization is negligibly small. The dashed line and dotted line in Fig. |21 d) show 
the split of ?/;(3686) peak due to the EM vacuum polarization given in Eq. (jHl). The label 
'bare' in Fig. |21d) indicate the lowest order peak position which can be seen in the processes 
in which there are no EM effects infected. 

The vacuum polarization effects also cause the flat continuum QED contribution vibration 
under the sharp resonance, which is shown by dashed dotted line in Fig. El d) and also in 
Fig. El b) by the red line The solid line in Fig. Eld) represents the expected total cross section 
in the energy region. 
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Figure 3: The simulation of a fine cross section scan in ?/'(3686) resonance region where the 
energy spread is assumed to be 1.35 MeV. In a), the expected observation cross section is 
shown by the sohd hne, while the green line shows the lowest cross section; and in b) the 
contributions from the three observation cross sections, i.e. the J/^ tail, the '0(3686) and 
the continuum QED process, are demonstrated. The circle marks on the lines in b) represent 
the integrated results of Eq. (jT)) for the cross sections of the three processes. 



Fig. El a) shows the different behavior of the lowest order cross section and the expected 
observable cross section due to radiative effects and with finite collision energy spread (1.35 
MeV) generated with the generator in the '?/;(3686) resonance region. While Fig. 01 b) shows 
the contributions of the observation cross sections of three relative processes in the energy 
region, the tail of J /ip., the '?/'(3686) resonance and the continuum QED, in which the vibration 
of the continuum QED line shape caused by the vacuum polarization corrections is still seen 
clearly in the finite energy spread case (the red line). The finite energy spread does drop the 
sharp resonance peak significantly but does not eliminate the vibration of continuum QED 
line shape as shown in the figure. If the data handling, the event selection, the determination 
of the detection efficiencies and the physical considerations of the hadron spectrum structure 
are all treated properly and correctly in the cross section scan experiments, one can expect 
that the vibration effects of continuum QED line shape due to vacuum polarization lying on 
the narrow resonance region would be observable. 

The circle marks on the lines in Fig. El b) are the integrated numerical results of Eq. (^J) 
for the corresponding cross sections. The accuracies of the simulated cross sections by the 
generator are quite good. 

2.2 Sub- generators 

From Fig. ^ one can see that the hadronic events not only be produced at the i/i, but also 
produced in the full energy range from y/s to 0.28 GeV, the two pion production threshold. 
To correctly simulate the full process produced in e^e~ annihilation, the generator should be 
a mixture of sub-generators in which each single generator simulates a special processes and 
runs together with others, driving according to the ISR return energy spectrum weighted 
by the lowest order cross sections of the processes involved in the energy region. Up to the 
charmonium system energy region, these physics processes include one photon continuum 
hadron production, e^e~ —>■ J/ f/', '0(3686), which decay to hadronic final states inclusively, 
the e~^e~ — ^ ^(3770), ?/'(4040) and '?/'(4160) etc. which go to the hadron final states through 
the mediate charmed meson pairs inclusively. The production of 1 resonances at lower 
energy range are also included, such as the p(1700), 0(1680), a;(1650), p(1450), Ci;(1420), 
0(1020), u;(780) and p(770) etc.. Those states are arranged respectively to go to various 
final states such as tttt, KK,KKn, A7r,57r,67r.... exclusively with their corresponding cross 
sections and branching fractions. The events of e~^e~ — ^ 7* — > hadrons and the inclusive 
hadron decays of charmonium resonance -0(3686) and J/ip are generated by their respective 
inclusive LUND-type generators jH] [11 • While the events which decay from other resonances 
such as the ^(3770), 0(1680), 0(1020), p(1700), p(1450), p(770) and u;(782) etc. to all 
possible final states according to the known decay modes and branching fractions j3] , can be 
generated by corresponding sub-generators conveniently. The Lorentz boost of the hadron 
system recoiling against the radiative photons are planted into the sub-generator's operation. 
A normal distribution random number generator has been employed to simulate the effects 
of energy spread at the collision energy ^/s. 



2.3 Sampling Technique 

From Eq. ((Tj), Eq. (0) and Fig. ^ one can see that there are many sharp peaks, such as 
the J/f/', ^^(SGSG) peaks (see Fig. [T]and Fig. |2)) and the sharp photon ISR spectrum at the 
singularity when x — *> (equivalent to y/s' — * 3.8 GeV in Fig. ^, which enter into the 
probability distribution function given in Eq. (^Q). All of those sharp structures may slow 
greatly the Mont Carlo sampling procedure and affect the accuracy of the event generations, 
especially for the accuracies of the shapes of the sharp variation lines which are extremely 
important for the determination of the detection efficiencies. Sometimes people can elim- 
inate the singularities by Monte Carlo inverse transformation techniques, but it is almost 
impossible to eliminate all those singularities simultaneously by an inverse transformation 
function. A special importance sampling method with multi-section scheme treatment has 
been developed in this work. One divides the whole of ISR return energy region into several 
sections to guarantee that for every section there is only one or less peak of the distribution 
function given by Eq. (^Q). Three kinds of sections are defined: 1) the resonance section in 
which there is only one sharp resonance like J /ip or '?/;(3686); 2) the slow varying sections in 
which there is no sharp resonance or any singularity; and 3) the section in which the ISR 
spectrum has the singularity i.e. the region which directly touches to the collision energy 
point X — s- 0, (s' — > s), for which the section range is variable and depending on the position 
of the nearest sharp resonance beneath the collision energy ^/s. When the collision energy 
with its energy spread range close to the narrow resonance or overlap each other, the width 
of the ISR singularity section xq (xq > x > 0) should ensure that both of the variation of the 
lowest order cross section in the singularity section region (xq > a; > 0) and the variation of 
the ISR spectrum in the touched sharp resonance section should be comparably moderate 
enough. 

The sharp resonance peak in the resonance section can be eliminated easily by the inverse 
function transformation, defining the new variant p instead of x, in Eq. (^ 

X = 1 - [M ■ r ■ tan{tan-\ ^ ^ ^ -n-p) + M^s, (12) 

M-V M-T ^ ' 

serves as the probability normalization constant in the transformation, and x„p and xiow are 
the up and low limits corresponding to the resonance section. The p is then the normalized 
random number taken uniformly from to 1, which keeps the event production probability, 
i.e. the integrated area of Fig.Qunchanged in the resonace section under the transformation. 
This transformation keeps the exact shape of the sharp resonance, but it blocks generating 
the initial state photons individually. It can just output the total energy and momentum of 
the emitted photons to serve for the Lorentz boost of the recoiling hadronic system. For the 
purpose of inclusive hadron measurement in e^e~ collision experiment this simplification is 
feasible. The price paid here is worthy for the precise measurement of the narrow resonance 
parameters. 

In this stage one can simply define the new normalized sampling random number ^ in 
(0,1) section instead of x to eliminate the singularity of the exponential ISR spectrum, the 



in which 



leading logarithm term of Eq. Q and keep the event production probability unchanged, 

X = x^e'^ (14) 

where the (5 is as defined in Eq. (0), and the section width xq here saves as the probability 
normalization constant in the transformation. 

After performing the importance sampling procedure for the multi-section process ac- 
cording to the distribution of Eq. ^ and the normalized ratios of the lowest order cross 
sections of the processes involved to drive the corresponding sub-generators the Mont Carlo 
events will be then generated finally as mentioned in above section. 
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Figure 4: The properties of hadronic events produced at a/s = 3.773 GeV. Shadows are the 
distribution of experimental data [H] , and lines are those generated with this generator; a) 
the charged multiplicity; b) the KNO scaling jTU]; c) thrust [TT], d) sphericity [12], e) apla- 
narity ^2j f ) Feynman scaling variable, g) pseudo-rapidity, h) and i) transverse momentum 
projection in and out the jet injection plane 
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Figure 5: The efficiencies for detection of tlie inclusive liadronic events at eacli nominal 
center-of-mass energies. 



3 Comparison with experiment data 

Properties of liadronic events can be described by distributions of their kinematic variables. 
Fig. 0] show the distributions of kinematic variables of hadronic events produced at ^/s = 
3.773 GeV. In the figure, shadows show the distributions of the variables from data jH], and 
lines represent the distributions from Monte Carlo sample generated with this generator, in 
which the detector simulation is based on GEANT3 package [H] . In Fig. IH a) is the charged 
multiplicity Nch of events; b) is the KNO scaling distribution ^2]; c), d) and e) are the 
thrust distribution ^Jj , the sphericity distribution and aplanarity distribution of events ^2] ■, 
respectively; f) is the distribution of the Feynman scaling variable. Kinematic properties of 
hadronic events are described by g) pseudo-rapidity , and transverse momentum distributions 
h), n) and i). From Fig. |3]one can see that the qualities of the general simulation to the 
data are satisfied. 



4 Application of the generator 

Feeding the reconstructed Monte Carlo events for e^e~ -^ hadrons generated at each of 
the nominal center-of-mass energy points, at which the cross section scans were performed, 
into the analysis program to select the inclusive hadronic events one can determine the 
detection efficiencies ehad{s). Fig. El shows the Monte Carlo efficiencies for detection of the 
hadronic events produced at the different nominal center-of-mass energies from 3.650 to 
3.872 GeV. From the figure one can see that the efficiencies at the different energies with 
different physics processes involved are quite different. The efficiencies determined with the 
generator were used in the measurement of R values around 3.650 and at 3.773 GeV [T^ 
which were used to extract the inclusive branching fractions of ?/'(3770) decay to DD, and in 
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the cross section scan experiment in the energy region covering both oiip{3686) and ^/^(STTO) 
resonances performed with BES-II detector at the BEPC coUider in Beijing ^3]. For more 
precise measurements of the resonance parameters obtained by fitting the cross sections for 
inclusive hadron production in the energy range, where the comphcate physical structure 
like the ip{3686), 'z/^(3770) and ip{A04:0) resonances involved, the correct determination of the 
detection efficiencies is crucial. 



5 Conclusion 

A new generator with a^ order radiative correction for the full simulation of e~^e~ -^ hadron 
process in the cross section scan experiment in the energy region covering multiple reso- 
nance structure is present. The complex resonances structure such as the ip{3770), '0(3686), 
J/ip, V'(4040) and ■?/'(4160) as well as the continuum QED production and the resonances 
with masses below 2 GeV/c^, are all included in the generator. The generator reproduces 
the properties of hadronic events very well and the detector efficiency can be determined 
conveniently for the hadronic cross section scan experiment. 
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